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Abstract: There is a deep cosmological mystery: although dependent on very differ- 
ent underlying physics, the timescales of structure formation, of galaxy cooling (both 
radiatively and against the CMB), and of vacuum domination do not differ by many 
orders of magnitude, but are all comparable to the present age of the universe. By 
scanning four landscape parameters simultaneously, we show that this quadruple co- 
incidence is resolved. We assume only that the statistical distribution of parameter 
values in the multiverse grows towards certain catastrophic boundaries we identify, 
across which there are drastic regime changes. We find order-of-magnitude predic- 
tions for the cosmological constant, the primordial density contrast, the temperature 
at matter-radiation equality, the typical galaxy mass, and the age of the universe, in 
terms of the fine structure constant and the electron, proton and Planck masses. Our 
approach permits a systematic evaluation of measure proposals; with the causal patch 
measure, we find no runaway of the primordial density contrast and the cosmological 
constant to large values. 
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1. Introduction 



In the evolution of the universe there are certain critical time scales that separate 
one era from another. The electroweak phase transition occurred at a time of order 
1CT 12 sec, the QCD phase transition at order 1CT 4 sec. and the transition from radiation 
to matter domination at order 10 12 sec. The time scale of inflation is very uncertain, 
and could be as long as a second or as short as 1CT 38 sec. Apparently the universe 
evolved smoothly for many decades in any given era, but between these eras, at well 
separated time scales, there were dramatic changes of regime. 

This picture breaks down in a remarkable way during the present era. There are 
several time scales that originate from very different fundamental physics, and would 
be expected to differ by many orders of magnitude, but all have the same order of 
magnitude. Each of these time scales by themselves could herald a change of regime. 
The change from matter to vacuum domination, instigating a new era of accelerated 
expansion, occurs at t\ A -1 / 2 , where A = 8-kG^Pa is the cosmological constant. 
Another change of regime occurs at the time t vir which is determined by a combination 
of the size of the primordial density perturbations and the temperature of matter 
radiation equality. At this time, density perturbations become non-linear. The universe 
transits out of the nearly homogeneous state it had for a number of eras and becomes 
clumpy. There are two possibilities for subsequent evolution. One is that the clumps or 
proto-galaxies cool very quickly to form galaxies, t coo i <C t vir , followed by fragmentation 
and star formation. In this case the change at t vir is from a regime of approximate 
homogeneity to one of star formation. The second possibility is that t coo \ ^> t vir , so that 
there is first an era of virialized proto-galaxies, and then a much later era during which 
they cool. In fact, two very different processes can contribute to cooling: radiative 
emission and inverse Compton scattering from the background radiation, with time 
scales t ra d and t CO mp, respectively, and the cooling time scale t coo i is given by the smaller 
of these two. Finally, once proto-galaxies have cooled, they give rise to additional time 
scales, including the time scale for the subsequent appearance of observers, t ODS , and 
the stellar burning time scale tbum- 

Other critical time scales in the past history of our universe span more than 60 
orders of magnitude. Those of future phenomena, such as the decay of matter, black 
hole evaporation, etc., can range over thousands of orders of magnitude [1]. One might 
expect that the five time scales £a, t vir , t ra d, t comp and t obs , dependent as they are 
on entirely different physics, would also be well separated. Remarkably, they are all 
comparable, 

t\ ~ tvir ~ ^rad ~ tcomp ~ tobs, (1-1) 
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having an order of magnitude of 10 10 years. 1 The picture of well separated eras breaks 
down during the present time, since many changes in regime are occurring simultane- 
ously. It is hard to accept that Eq. is purely coincidental; but it is also hard to 
see how it could be derived from symmetries of some underlying field theory. 

If A varies in the multiverse, one of the coincidences of Eq. (|1.1|) can be explained 
by environmental selection. Only those universes with t& > £ v ir have any large scale 
structure and can be observed [2] (see also [3]). Given that most universes have a large 
cosmo logical constant, the statistical prediction over the multiverse is t& ~ t v[r . This re- 
sult seems particularly important since it is the only known solution to the cosmological 
constant problem, and it predicted the discovery of nonzero vacuum energy. In fact the 
observed value of the cosmological constant is somewhat smaller than expected from 
this result [4], but it is important to note that the prediction depends on the measure 
used to define multiverse probability distributions. Using the causal patch measure [5], 
the prediction is modified to t& ~ t v ; r + t coo i + t b s , and is more successful [6]. 

It is important to understand that it is legitimate in such calculations to hold all pa- 
rameters but one fixed, even though many parameters can vary in the string landscape. 
After all, it yielded a first, nontrivial prediction that could have failed (e.g., if A had 
not been discovered). Once additional parameters are allowed to vary simultaneously, 
the theory has another opportunity to fail, or to succeed. A long-standing concern 
about the prediction for the cosmological constant is that simultaneous scanning of 
the virial density leads to runaway behavior favoring large values of both quantities. 
However, runaway behavior need not ensue if the range of new scanning parameters 
is constrained by additional catastrophic boundaries, and if multiverse forces point 
towards these boundaries. 

In this paper, we study this possibility. We identify a number of possible catas- 
trophic boundaries. We argue that under suitable assumptions about the prior distri- 
bution of parameters in the landscape, these boundaries can prevent runaway. What 
makes this viewpoint especially compelling is that it explains all of the coincidences in 
Eq. fll.lQ , allowing for novel cosmological predictions from environmental selection. 

Outline In section |2|, we review how predictions arise from multiverse probability dis- 
tributions that grow towards catastrophic boundaries, across which there are significant 
changes in the anthropic weighting factor [7]. We stress the utility of the multiverse 
force, the logarithmic derivative of the distribution function, and of a set of variables 
chosen to be orthogonal to the catastrophic boundaries. In section |3] we review analytic 

*As a by-product of understanding these coincidences, we will find that the time scale for a main 
sequence star of maximal mass, ibum, is constrained to be parametrically equal to these other time 
scales. 
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estimates for a number of astrophysical processes: halo virialization, halo disruption 
by a cosmological constant, radiative cooling, and inverse Compton cooling of proto- 
galaxies. Our aim is to obtain an approximate parametric dependence of astrophysical 
scales in our universe on parameters such as the fine structure constant, a, the electron 
mass, m e , the proton mass, m p , and A, which may vary in the multiverse. 

In sections f|, |5], ^ and [F] we successively allow the virial density, the cosmological 
constant, the time scale of observers, and the Compton cooling time scale to scan. In 
each section, we add one more scanning parameter and one more catastrophic boundary, 
and we obtain one new prediction, while maintaining the previous predictions and 
results. At each stage, we derive the range in which multiverse forces would have to lie 
so that our universe is not atypical. 

In section [|, we assume that the matter density at virialization scans in the multi- 
verse. We identify a catastrophic boundary: for sufficiently small density, proto-galaxies 
cannot cool. Assuming a multiverse force towards small density, we find predictions for 
the virial density and the mass of a galaxy, 



a 4 



a 5 

M S-l ~ n2 1/2 5/2 > ( L3 ) 

G^m e nip 

in terms of a, m e , m p , and the Newton constant, G^. (We have omitted numerical 
coefficients that are displayed in the main text.) This also explains the coincidence 
that the timescales for radiative cooling and galactic halo virialization are comparable: 

trad ~ t vil . (1.4) 

The distance of the observed virial density from the catastrophic boundary is used to 
place empirical bounds on the strength of the multiverse force. 

Next, in section we consider scanning the cosmological constant A in addition 
to the virial density. In a certain class of measures that includes the scale factor cutoff 
measure proposed in Ref. [8], A is constrained from above by the disruption of halo 
formation. There is a potential runaway favoring large values of both quantities, since a 
larger virial density allows for a larger cosmological constant while preserving structure 
formation. However, for a certain range of multiverse force directions, this runaway 
is avoided and one predicts that typical values of the two parameters are close to the 
intersection of the two boundaries. This leads to a second prediction, now for the 
vacuum energy: 

1 G^rr^m 2 

^ ~ TTTT — 4 • ( L5 ) 

Gn^a or 
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It explains, moreover, the double coincidence 



t\ ~ ^rad ~ t v ir- 



(1.6) 



As in the previous section, we derive empirical constraints on the values of multiverse 
force components by comparing the observed pa and p vir to the critical values. We can 
compare these constraints to theoretical predictions, which depend on the measure on 
the multiverse, and on the prior distribution in the landscape (which is known for p\). 
For the measure of [8], the empirical and theoretical limits conflict. This illustrates 
how our approach can be used as a systematic discriminator between measures. 

In section |^, we consider a different class of measures, where the weight of observa- 
tions made after t\ is exponentially suppressed. (This class includes the causal patch 
measure [5], the causal diamond measure [6], and the modified scale factor measure 
of [9].) In the sharp boundary approximation, catastrophe occurs unless t\ > ^vir + ^obs- 
This condition is approximately equivalent to two boundaries — the upper bound on 
the cosmological constant that appeared in the previous section, and a new boundary 
^obs < ^a- Correspondingly, this class of measures allows us to scan a third parameter, 
tobs, in addition to scanning p v i r and pa- A wide range of multiverse forces will favor 
universes close to the intersection of the three catastrophic boundaries, allowing an 
explanation of the triple coincidence 



We determine quantitative empirical constraints on force components and find that they 
are perfectly compatible with the prior distribution of pa if the causal patch measure 
is used. 

Halos that virialize in a sufficiently hot universe cool by inverse Compton scattering 
against the cosmic background radiation. We do not know whether this change of 
regime is catastrophic, but in section [7], we explore the consequences of assuming that 
it is. With this new boundary, we can allow a fourth parameter to scan, which we 
choose to be the strength of primordial density perturbations Q. In this case a further 
prediction follows from assuming a multiverse distribution that prefers large values of 



Since Q and p v ; r determine the temperature at equality, combining this with Eq. ( |1.2|) 
yields 



iobs ~ ^A ~ irad ~ ^vir- 



(1.7) 



Q 






a 
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This additional prediction, combined with Eq. (|1.7|) , explains the quadruple coincidence 

Wip ~ ^obs ~ t A ~ t rad ~ £ vir . (1-10) 

In section |8], we discuss how far this program of adding more scanning parameters 
and catastrophic boundaries can be taken. We argue that many more predictions are 
possible, and we consider various ultimate solutions to runaway behavior. 

2. A Multiverse Force and Unnatural Predictions 

In this section, we review the approach to multiverse predictions proposed in Ref. [7]. 
We will also discuss the quantitative relation between the strength of a multiverse force 
and the proximity of typical observers to the corresponding catastrophic boundary. 

A landscape of vacua, such as the string landscape, will contain scanning param- 
eters, x, that vary between vacua [10]. We can consider any parameters that enter 
physical phenomena; they need not be Lagrangian parameters. The multiverse distri- 
bution 

dN = f(x)d\nx (2.1) 

describes the relative number of times the value x is observed, and thus its probability. 
In most situations, the distribution 

f{x) = f(x)w(x) (2.2) 

factorizes into an a priori distribution, f(x), of vacua in the theory landscape and an 
anthropic weighting factor, w(x), roughly the number of observers in vacua where the 
parameter takes on the value x. (This number naively vanishes or diverges for all x 
because of eternal inflation, and a cutoff or measure is required to make it well-defined.) 

Living on the edge In general, to make a prediction one needs to understand both 
f(x) and w(x) — the landscape, the measure, and the abundance of observers — in great 
detail. Under the following two conditions, however, quantitative predictions can be 
made even though only qualitative features of / and w are known: 

• There are catastrophic boundaries where w(x) drops sharply. Studying physics 
in the parameter space of x, one often discovers that there are boundaries in this 
space across which there is a sudden change in the physics. For example, this 
could be a phase transition or the stability criterion for some complex object. In 
the extreme limit, w(x) can be represented by a theta function, so that observers 
are found in universes on one side of the boundary, but not in universes on the 
other side. This limit need not be exactly realized in practice, but it will simplify 
our analysis considerably. 
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x b 

Figure 1: A strong multiverse force towards a catastrophic boundary at cc& implies that a 
typical observer will be located close to the boundary at x a . 



• In the region of x where observers are possible, the probability distribution, f(x), 
increases towards the catastrophic boundary(s). Thus, most universes do not 
support observers. 

In this situation, a typical observer would expect to measure x close to the boundary. 
Vacua with very different values of x either contain no observers, or are very rare in 
the multiverse. 

The crucial point is that this quantitative prediction depends mostly on the location 
of the boundary, which can be calculated using the conventional law of physics. It 
depends only weakly on the multiverse, in that only the above qualitative assumptions 
on the multiverse distribution must be satisfied. This situation is illustrated in Fig. |l| 
for the case of a single scanning parameter with a catastrophic boundary at x = Xf,. 

The prediction x ~ x^ does not follow from any symmetry of the theory; indeed 
it may appear coincidental. For example, it may relate two time scales or two mass 
scales that do not have a common origin. Since there is no symmetry reason that the 
observed value x Q should be close to Xb, from the conventional viewpoint the multiverse 
predictions appear quite unnatural. 

Taking x to be the vacuum energy density, p\ = A/87tGn, Fig. [I] illustrates 
the prediction of the cosmological constant, with Xb ~ pvn- Using a logarithmic 
scale, the multiverse distribution indeed grows towards the boundary. This repre- 
sents the cosmological constant problem. The naturalness probability, defined by 
P = | f* b f(x)dx/ ff(x)dx\, is of order 10 -120 , and from the usual symmetry view- 
point, it is extremely unnatural for p\ to take the observed value. In the multiverse 
picture, however, the vast majority of universes are simply irrelevant since they contain 
no large scale structure and thus no observers. Observing p^ ~ p vir , therefore, is not a 
mystery. 

Force strength The steepness of the distribution function f(x) determines how close 
a typical observer lies to the boundary. We can refine our analysis by defining the 
multiverse probability force 

din/ , s 

F=— '-. 2.3 
a In x 
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In this paper we will consider only power law distributions, f(x) oc x p . Then the force 
is simply given by the power, F = p. Defining x v to be the point such that the fraction 
of observers farther from the boundary than x u is u, we find 

^ = z//p. (2.4) 

x b 

Thus x . 5 /xb is a useful measure of the typical "distance" of an observer from the 
boundary; it is 1/2 for p = 1 and rapidly becomes smaller (larger) as p decreases 
(increases), since 1/p appears in the exponent in Eq. (|2.4j) . For p = the distribution 
is logarithmic, with equal population of all logarithmic intervals. In this case, one does 
not expect to find x close to x b , so it is of no interest here. 

If the force is somehow known, then we can predict not only that x should be close 
to Xb, but how close it should be. Roughly, one expects to measure a value of order 
X0.5, and one would be mildly surprised to measure values farther from the boundary 
than, say, £o.i6, or closer than xo.84- 

In practice, the strength of the force is not precisely known from the theory of the 
landscape. We can, however, still argue that the multiverse is likely playing a role in 
determining the value of a parameter x if x Q is anomalously close to x^, since such a 
phenomenon is extremely hard to explain in any other framework. In this circumstance, 
one can reverse the logic given above and derive bounds on the strength of the force; 
for example, if x is extremely close to Xb, the force is likely to be quite strong. This 
information can then be used to study if a certain assumption on the landscape or 
measure is consistent with data. 

Let us make this point more precise. Given x Q and x b , the observed and boundary 
values of x, Eq. ( |2.4[ ) allows us to compute the force p such that a fraction v of observers 
are farther from the boundary than we are: 

= £&■ (2 ' 5) 

For example, suppose that we are willing to reject models that would put us farther 
than one standard deviation from the median. Then we can infer that the multiverse 
force lies in the range p(0.16) to p(0.84). If we allow two standard deviations, the force 
must lie in the range p(0.025) to p(0.975). 

As we shall see below, such empirical bounds on the multiverse force can be quite 
useful, for two reasons. First, there may be more than one way of obtaining a bound 
on the multiverse force acting on a particular parameter, allowing for a nontrivial 
consistency check. Second, even though we may not be able to compute the multiverse 
force from first principles, there may be theoretical arguments why it should satisfy a 
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X2 A 




Figure 2: With two scanning parameters, (xi,x%), a multiverse force towards a single catas- 
trophic boundary, such as A, may lead to runaway behavior. The introduction of a second 
boundary, B, can prevent this runaway and lead to predictions for both parameters. 

certain upper or lower bound. In these cases, comparison of the theoretical arguments 
to the empirical bounds offers a possibility of falsifying various assumptions on the 
properties of the landscape and the measure. 

Multiple parameters Suppose that two parameters scan. Consider, for example, 
the catastrophic boundary labeled A in Fig. ^|. In general the multiverse force is not 
perpendicular to the boundary, so there is now a runaway problem, as shown by the 
dashed arrow. This runaway behavior may be halted in several ways. A special shape 
of the boundary could prevent runaway, but for many boundaries this is not the case. 
The runaway could be prevented by f(x) reaching a maximum at some point along 
the boundary; however, in this case this special point on the boundary results from an 
assumption on the form of f(x) and is not predicted. 

Alternatively, the runaway problem can be solved by introducing a second catas- 
trophic boundary, as illustrated by the boundary B in Fig. ^. In this case, the bound- 
aries define a special point by their intersection, and a multiverse force in a wide range 
of directions will lead to typical observers being close to this point, as shown by the dot 
in Fig. ||. Not only does the second boundary stop the runaway behavior, it leads to a 
second prediction, which will again appear coincidental and unnatural from a symme- 
try viewpoint. Typical observers are in universes near the "tip of the cone", a special 
point in two-dimensional parameter space. 

With n scanning parameters there could be many runaway directions. But with 
n boundaries of codimension one, it is possible that all the scanning parameters are 
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determined. With the multiverse force pointing in a generic subset of directions, the 
most probable observed universes will be at the tip of an n-dimensional cone, giving n 
predictions. 

In general, these predictions will relate the scanning parameters x to fixed parame- 
ters y, which are not allowed to scan — either because they do not vary in the landscape, 
or because we choose to consider only a subset of the landscape. It is possible that, 
for the phenomena of interest, all parameters scan, and the configuration of bound- 
aries together with our assumptions on f(x) prevent any runaway and determine all 
the underlying parameters. In practice it may be hard to come up with such a closed, 
completely determined set of parameters, and the multiverse may not lead to such sit- 
uations. Still, it is worth exploring how far this program of multiverse predictions can 
be taken. 

In this paper, we consider mainly power law forces, f(xi,..., x n ) oc Y%i=i $t ■> an d 
power law catastrophic boundaries, YYj=i ^ K = A? * = which generically 

intersect at a critical point (xi jC , . . . , x njC ). It will be convenient to work with the 
logarithmic quantities 

Zi = \n^, (2.6) 

and to use capital letters Z, P, and B, to denote the vectors and matrix constructed 
from Zi, Pi, and bij. Then the critical point is at Z = 0. 
Let us define new variables Ui = b^Zj, or 

U = BZ, (2.7) 

which measure the distances orthogonal to the n boundaries. The boundaries are 
defined by U{ = 0, i = l,...,n. This is illustrated in the two-dimensional case in 
Fig. Without loss of generality, we will choose signs such that Ui > in the allowed 
region. The variables U are convenient to work with because the power-law form of 
the distribution function is preserved if one or more of the Ui are integrated out. This 
is not generally the case for the X{ or Z{. 

With the new notation and variables, the distribution function takes the form 

\nf(Z) = P-Z = P-B- 1 U = S-U, (2.8) 

where B~ x denotes the inverse matrix, and we have defined 

S = (B- l ) T P. (2.9) 

The force vector, with components pi = in the original coordinates, has compo- 
nents Si = orthogonal to the boundaries. Importantly, Sj is the effective one- 
dimensional force on the variable Ui if all other Uj, j ^ i, are integrated out. (Note 



-10- 




ln^i 



Figure 3: Given n boundary hyperplanes, it is convenient to work with new variables ut 
such that the allowed region is a quadrant and the critical point is at the origin of R n . In 
these variables, the distance from each boundary Ui is inversely related to the corresponding 
multiverse force component Sj, with no mixing between components. 



that the analogous statement would not hold true for the original forces and variables 
Pi and Ui.) 

To see this, consider the one-dimensional distribution functions fi(ui) = which 
can be derived from the multivariate distribution by integrating out all but one variable: 



POO 

fi(ui) oc / Y\ duj I det B]" 1 exp(5' ■ U) oc exp(sj-Uj). 



(2.10) 



The one-dimensional effective force on m is 



— = Si, (2.11) 

dUi 

as claimed. With our convention that the Ui are positive in the allowed region, the 
must all be negative to avoid runaway. 

The fraction of observers farther than Ui from the boundary m = is 

^ = exp(siUi). (2.12) 

If the multiverse force is known, the characteristic distances Wj can be predicted in 
analogy to Eq. fl2.4Q . For example, from Ui = 0.5 one obtains the median values 
Ui = s" 1 In 0.5. With Z = B^U and Eq. (|2.6f) , this translates into predictions for 
characteristic values of the original parameters Xj. 
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Conversely, if the multiverse force is not known, it can be constrained by demanding 
that the observed values not be improbably close to or far from the boundary, as in 



These n constraints can be translated into constraints on the n original force compo- 
nents pi, with P = B T S. 

The region of parameter space that allows observers - the observer region - is 
assumed to be open in the above analysis. However, there could be more catastrophic 
boundaries that make the observer region finite. If these additional boundaries are not 
distant from the point describing our universe, the analysis would change accordingly. 

3. Parametric Dependence of Astrophysical Scales 

In this section, we will discuss astrophysics in our universe. We will ask how astro- 
physical quantities, like the virial density and temperature, and the range of galactic 
masses, are related to "fundamental" parameters such a, m e , m p , Q, etc. Such calcula- 
tions have nothing to do with the landscape or with imposing anthropic cuts; they are 
merely derivations of observed macroscopic phenomena from observed microphysics. A 
large enough computer could perform them with arbitrary accuracy, given the observed 
values of the input parameters. 

We will not attempt a very detailed analysis. Since we will be mainly interested in 
the parametric dependence on the fundamental parameters, our goal is to capture the 
dominant physics that controls the phenomena we study. In order to compensate for 
analytical shortcuts and bring our simple models into precise agreement with data, we 
will insert suitable order one factors where necessary (for example, c lon and c vir below). 

Note that this does not amount to putting in the answer to the central problem 
investigated in this paper. Our goal is to explain why certain combinations of funda- 
mental parameters (namely those that appear in derived astrophysical and cosmological 
quantities) take on the values they do. This question becomes meaningful only in the 
larger setting of the multiverse, where what would otherwise be fixed input parameters 
like a, m e , m p , Q, etc., can vary. It will be considered in the following sections. 

3.1 Halo formation 

Complex structures, such as stars, will not form unless the virialized gas inside dark 
matter halos is able to cool. In our universe, there is both a maximum and a minimum 
characteristic halo mass for which radiative cooling is successful. (Compton cooling 
does not play an important role in our universe and will be considered in section |7|.) In 



Eq. @: 




(2.13) 
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this section we review the conditions for cooling and derive the upper and lower bounds 
on the halo mass. 

We stress that the present section is entirely about reviewing the conventional 
physics of cooling in our universe; "fundamental" parameters like Q or m e are taken 
to have the observed values. In later sections, we will consider multiverse settings in 
which one or more parameters scan. At that point we will make use of the parametric 
dependence of the expressions we derive here, in order to explain why certain parameters 
that scan over the multiverse are observed to take the values that they do. 

We assume adiabatic, near scale-invariant, density perturbations with primordial 
amplitude Q. Far into the matter dominated era, a density perturbation of mass M 
grows linearly with the scale factor, a(t), and at time t has an amplitude 

Q(M,t) = (rt±Y 3 f(M)Q, (3.1) 



Peq 



where p(t) is the matter density and 



«- - (3-2) 

is its value at the time t eq when it equals the radiation density. Here, N eq ~ 3.36 is the 
number of relativistic species at the temperature of matter-radiation equality, and T eq 
is the temperature at equality. We will consider T eq the fundamental parameter that 
defines the time of equality. 

The function f(M) includes both a mild logarithmic growth for perturbations that 
entered the horizon before t eq and complicated behavior over a lengthy transition region 
between radiation and matter dominated eras. For example, for M = 1O 12 M , which 
is close to the mass of the Milky Way, f(M) ~ 43 [11]. 

For sufficiently small p^, the perturbation on scale M goes non-linear and virializes 
at time defined by 

Q(M,t vir ) = 5 col = 1 2Q ; ~ 1.69. (3.3) 
With Eq. flOP, one finds 



tvir = Cvir U^cJW g7p) ' (3 - 4) 

where the coefficient c va ~ 0.5 is introduced to account for the effects of halo merging. 2 
We have defined a fiducial density 

P = Q 3 T e % (3-5) 



2 Halos grow by mergers and by accretion. Both effects contribute to the mass M in Eq. (|3.3[), and 
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which is closely related to the virial density of the halo 

n f_P ir 2 N e J p f(Mf 3 4 _ Tr 2 N eq f p f{M) 3 

Pvir " 6vrG N e r " 30 QI ^~ 30 <^<4 P ' ^ 

in an approximation where the virialized halo is taken to be spherically symmetric with 
uniform density. Note that p v ; r is a factor f p = 187r 2 larger than the ambient density 

cltj ^ vir • 

3.2 Radiative cooling and galaxy formation 

Galaxies form in halos only if the virialized baryonic gas can cool and condense. The 
temperature of the virialized halo is set by the gravitational potential energy gained 
by electrons and protons falling into the halo. By the virial theorem, 

1 SGnMm = 3 (37) 

2 5.iLyir 2 

We will take the average molecular mass /i to be m p /2, so that 

5 V 3 / 10*coi V«^,J V ' 

In our universe, Q ~ 2.0 x 10 -5 and T eq ~ 0.82 eV, so that a halo of mass 1O 12 M has 
t vir ~ 3.6 x 10 9 years, p vir ~ (1.5 x 10~ 2 eV) 4 and T vir ~ 40 eV ~ 4.6 x 10 5 K. 

We can now state the conditions for cooling. First, the temperature of the halo 
must be large enough to give significant ionization 

T vir > c ion a 2 m e , (3.9) 



e — 



where, for bremsstrahlung and hydrogen line cooling, C; on ~ 0.05 so that 
1.5 x 10 4 K. Using Eq. (|3.8|) , this ionization requirement leads to a minimum value of 
the halo mass 

~ V nf(M) ) {2N cqfp ) GfrnTpV* ^ 

We will comment at the end of this section on the numerical value obtained here. 



they cannot be disentangled in a simple spherical collapse model. However, only major mergers lead 
to virialization. In a more careful analysis, the density and temperature of the virialized gas is set not 
at the time when the halo reaches mass M, but at the earlier time when it underwent its last major 
merger. Our choice of c v i r mimics this effect; we thank S. Leichenauer for help with selecting the value 
0.5. 
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Second, the time scale for radiative cooling of the halo, t ra d, should be less than 
the dynamical time scale td yn [12-14], which we choose to be t vir : 

If T vir > 10 6 K, the halo is fully ionized, and the radiative cooling rate is dominated by 
bremsstrahlung emission with a time scale at virialization of 

n 3/2 rrl/2 

9 m p T v { T 

''brems „ T ^ o j l°- ±z J 

oK a A pB 

where K ~ 3.17, and ps = Jb Pvn is the baryon density in the virialized halo with 
Jb — 1/6. At lower temperatures other processes are important, for example hydrogen 
line cooling for 10 4 K < T < 10 6 K, so that we write the radiative cooling rate as 

n \ 3/2 rpl/2 



W = ( ^ / rad (T/m e , a) ) — 3 p vir . (3.13) 
\8K J ol pb 

The cooling condition of Eq. ( |3.11| ) is satisfied only by halos with mass 

M ~ M += 3^3 .2 f 3 ) - 3 9 /2 9/2 - ( 3 - 14 ) 

\ ° °col C virirad / G^m e m p ' 

We have thus translated the ionization and cooling rate conditions into lower and upper 
bounds on the halo mass, Eqs. ( |3.10| , p.!4| ). 



In our universe, M_ ~ 2.9 x 1O 9 M and M + ~ 4.3 x 1O U M . This interval is 
somewhat smaller than the observed range of halo masses; in particular, it does not 
include our own halo, with M ~ 1O 12 M . This has a number of reasons. We have 
worked with la overdensities to associate a typical virial density and temperature to a 
given mass scale. But a reasonable fraction of halos will form from, say, 2a overdensities. 
This can be incorporated into Eqs. ( EUOl EM) h Y substituting f(M) 2f(M); it 



lowers the minimum mass by about a factor of 3 and raises the maximum mass by a 
factor of 8. 

With this definition, the range does include the Milky Way, if only barely. This is 
as it should be: a halo of mass 1O 12 M is only just able to satisfy the cooling crite- 
ria [15], corresponding to the coincidence that t coo i ~ i vir . The lower bound, however, 
is still about two orders of magnitude higher than the smallest galaxies observed in 
our universe, which have M w 1O 7 M . This is because we have not included molecular 
cooling in our analysis. Molecular cooling can lower the cooling boundary by two orders 
of magnitude in our universe, but it will not have a significant effect in the neighbor- 
hood of the critical value of p which we will find in the following section. Hence we will 
ignore it here. 
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3.3 Disruption of halo formation by a positive cosmological constant 

For definiteness, we will consider only positive values of the cosmological constant in 
this paper. Negative values require a different analysis, though our ultimate conclusions 
would remain essentially unchanged. 

A cosmological constant, A > 0, or vacuum energy, p\ = A/87tGn, will come to 
dominate the evolution of the scale factor when the universe reaches an age of order 
— (3/A) 1 / 2 , when it begins to disrupt the formation of large scale structure. At late 
times, it causes the universe to expand exponentially with characteristic timescale t\. 

A primordial overdensity which, in a flat universe with vanishing vacuum energy, 
would evolve into a halo with virial density p v ; r will fail to evolve into a gravitationally 
bound object in a flat universe with [2] 

pa > |s. (3.i5) 

Using Eq. (|3.6|) , one finds that the observed vacuum energy, 

Pa,o = (2.3 x KT 3 eV) 4 (3.16) 
prevents the formation of halos of mass greater than 5.4 x 1O 14 M in our universe. 

3.4 Compton cooling 

A proto-galaxy can also cool by inverse Compton scattering of electrons from the cosmic 
background radiation. The cooling rate of a fully ionized plasma with temperature T 
against a background with temperature Tcmb <C T is given by 

t = ic^U^ 
1 Sm e 

where <tt = Qua 2 /?>m 2 e is the total Thompson cross section and a = 7r 2 /15. Note that 
this rate does not depend on the density or temperature of the plasma. Its inverse 
defines the Compton cooling time scale, 

135m 3 

^comp = TJTj q 2T"4 ' (3.18) 

which depends on cosmological time t through the CMB temperature 

T CMB (t) = 1 (3.19) 

7r 4 A eq G N i eq r 
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Compton cooling is effective as long as this time scale, at the time of virialization, 
t v i r , is faster than the dynamical time scale of the virialized halo, which is also given 
by t vir : 

^comp(^vir) ^ ^vir- (3.20) 

Thus, Compton cooling is effective for any halos that virialize prior to the time 

W„ - 39/57r7/5Ar 4/ 5 G 4/5 T 4/5 m 9/5 " ( 3 ' 21 ) 

In our universe, this time is approximately 0.36 Gyr. 

Even before this time, Compton cooling may not be the fastest cooling process. It 
will be faster than radiative cooling only if t com p < ^rad- However, we shall not need 
to investigate this condition separately in this paper, because we will be interested in 
situations where Compton cooling is effective for halos in which radiative cooling is on 
the verge of being ineffective, i.e., for halos with t vir ~ t rad . 

4. Predicting the Virial Density and Mass of Galactic Halos 

In this section, we will consider the consequences of allowing the characteristic virial 



density, p, to scan in the multiverse. In section |4.1| , we use the parametric expressions 
for the two cooling boundaries on the halo mass M derived in the previous section to 
show that there is a catastrophic lower bound on p, where the maximum and minimum 
mass become equal. We assume a multiverse force that pushes p towards the catas- 



trophic galaxy cooling boundary. In section |4.2| , we use this qualitative assumption to 
obtain rough quantitative predictions of the mass of a typical halo and its virial density, 
in terms of a, m e , m p , and G^. In section fO| , we compare our result to astrophysical 
data. As expected, our universe is near but not exactly at the catastrophic boundary. 
The distance to the boundary allows us to derive empirical bounds on the strength of 
the multiverse force. 

4.1 The catastrophic cooling boundary 



Both the upper and lower bound on halo masses, Eqs. Q3.10| , |3.14| ), depend on the 



primordial density contrast and the temperature at equality only through the combi- 



nation p = Q T e q, which can be thought of as the characteristic virial density in the 
universe. We will now assume that p is a variable that can take on different values in 
the multiverse. The origin of this scanning can be via Q or T eq , but we need not specify 
which in the analysis here. For now, we will regard Eqs. (|3.10| , p,14j ) as functions of the 



scanning parameter p, and we will phrase our assumptions in terms of the multiverse 
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Figure 4: (a) Simplified sketch of the boundaries for ionization (lower) and radiative cooling 
(upper) for a halo of mass M. Our galaxy, shown by the dot, is very close to the radiative 
cooling boundary and has a virial density only 1-2 orders of magnitude larger than the 
minimum possible, (b) Dotted lines show modifications to the ionization and radiative cooling 
boundaries from molecular cooling and atomic line cooling. To the right of the dashed line, 
halos are able to cool by inverse Compton scattering from the background radiation. 



distribution function for p. (Later, in section [7], we will consider an additional boundary 
that distinguishes between Q and T eq .) 

These functions are sketched in Fig. [|a (with additional cooling boundaries given 
in Fig. [|b). Apart from the mild logarithmic factor f{M), they are both simple power 
laws, one increasing and one decreasing with p. They cross at the critical value 

p c = K (0.97 x 10- eV)'. (4.1) 



It is important to note that M is not a scanning parameter. For any value of p, halos 
with all values of M will virialize. 3 We will assume, for simplicity, that observers form 
as long as there exists some range of M, no matter how small, that can cool. Thus, p c 
represents a catastrophic boundary in p space. Universes with p > p c contain observers; 
universes with p < p c do not. 

At the catastrophic boundary, only halos with critical mass 

M c = M + (p c ) = M.(Pc) - f 40 ^^; 1011 ^ ) n2 f /2 5/2 - 2.4 x 1O 1O M (4.2) 



3 This assumes that curvature and vacuum energy are negligible. We will consider nonzero vacuum 
energy in later sections. 
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and critical virial density 



V 



'vir.c - I 95 f K2 1 2 / ft4 



1-^ 



(7.7 x 1(T 3 eV) 4 (4.3) 



are able to cool. 

Although we have been keeping track of numerical factors in the above equations, 
we aim only to define the position of this catastrophic boundary at the order of magni- 
tude level. For example, there is considerable uncertainty in the numerical factor in the 
inequality of Eq. ( |3 . 1 1| ) . The physics of cooling involves non-uniformities and shocks 
and, furthermore, the boundary is not completely sharp, with the statistics of merging 
playing an important role. Moreover, as discussed in the previous section, there is some 
ambiguity in how M + and M_ should be defined. These ambiguities arise in part be- 
cause we are treating the catastrophic boundary as sharp; in a more detailed analysis, 
the anthropic weighting factor would change continuously over a range of values of p. 
Even then, a precise estimate of the number of observers would remain a challenge. 

4.2 Predictions 

To obtain a prediction, we assume that the multiverse distribution f(p) is a decreasing 
function. Thus, multiverse statistics favors low values of p, as shown by the direction 
of the arrow in Fig. [|a. In this case, observers typically measure a value of p close to 
the critical value: 

P ~ Pc- (4.4) 

Most universes have lower values of p, so that no perturbations are able to cool; these 
universes will not be observed. Observers live only in universes with p sufficiently large 
to allow cooling. Most will find themselves near the catastrophic boundary, with p 
larger, but not much larger, than p c given in Eq. Q4.3|) . 

There are a number of corollaries to this prediction. Galaxies in these universes 
will arise in a relatively narrow range of halo masses, controlled by the critical value 
M c defined in Eq. (FO), and ranging from a minimum halo mass of order M_ « 



M c (p/p c ) -1 / 2 to a maximum halo mass of order M + rs M c (p/p c ). (The minimum mass 
can be lower with molecular cooling, and we have neglected logarithmic factors.) In 
these halos, virial densities will be larger, but not much larger, than the critical value 
given in Eq. ( |4.3|) ; and the timescale for radiative galaxy cooling will not be much 
shorter than the timescale for virialization. 

For the sake of clarity, let us summarize these results by restating the relevant 
formulas without logarithmic or order one factors. By allowing Q and/or T eq to scan, 
and assuming a multiverse force towards small values of p = Q 3 T^ q , we were able to 
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predict the observed virial density at the order of magnitude level from fundamental 
parameters: 



a 4 



(4.5) 



It followed that our universe should be close to the radiative cooling boundary of 
Eq. (gXjj) : 



a 2 

irad ~ ^vir ~ ~ n , (4.6) 

explaining one of the coincidences of Eq. (|1.1|); and that the galactic mass scale arises 
from 4 

a 5 

M ~ (4 7) 

n 2 1/2 5/2- V*-^ 

G^rrie m p ' 

Another corollary is worth pointing out. The radiative cooling time scale at the critical 
point, Eq. (|3.11| ), has precisely the same parametric form as the lifetime for a main 
sequence hydrogen burning star of maximum mass, tbum- Hence, in a typical universe 
that is not far from the critical point, there is no freedom to choose the stellar lifetime: 

^burn ~ ^rad ~ ^vir- 

These predictions are successful. In our universe, p Q = QlT 4 q o ~ (2.4 x 10~ 4 eV) 4 , 
which is only a factor of 40 larger than the critical value: 

Po = 40 p c . (4.8) 

(We have restored order one factors.) This is remarkable since these two quantities could 
have been many orders of magnitude apart. It would appear even more impressive had 
we compared quantities of mass dimension one: pl^ 4 = 2.5 p l J A . Nevertheless, p Q is 
sufficiently larger than p c that in our universe galaxies with a considerable range of 
mass can cool radiatively. Our galaxy is near the upper end of this range, suggesting 
further environmental selection in our universe. 

4.3 The strength of the multiverse force 

The strength of the multiverse force determines how close a typical universe is to 
the critical point. For p, this strength cannot be predicted from first principles with 
current theoretical technology. Instead, let us now reverse course and put observational 
constraints on the force. From Eqs. ( |2.5| , f4.8|) , we conclude that the force satisfies 

-vp = o-i9 + -oit(^)- (4.9) 

4 Arguments that radiative cooling introduces a critical mass scale for galaxies, given by Eq. ([4.21). 
were made in the late 1970s [12-14, 16]. Since a large range of halos on one side of the boundary 
can cool, however, a multiverse force and anthropic arguments are needed to explain why observed 
galaxies are so close to this critical value. 
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Note that the force discussed here is the effective force obtained after integrating 
out parameters other than p [7], such as the cosmological constant. In the following 
sections we will consider such parameters, and we will be able to determine under which 
conditions a force in the above range can arise. We will see that this places interesting 
constraints on the multiverse distribution function, and in particular, on the measure 
for regulating the infinities of eternal inflation. 

5. Predicting the Cosmological Constant 

In this section, we will allow a second parameter to scan: the cosmological constant 
or vacuum energy, pa- 5 To avoid runaway, we need a second boundary. One way to 
obtain such a boundary is to require, as in the previous section, that galaxies form [2]. 
This condition is clearly necessary; whether it is the most important constraint on p^ 
depends on the choice of measure. As shown in Ref. [9], the version of the scale factor 
measure formulated in Ref. [8] leads to this constraint. 

We will be able to explain the coincidence t v ; r ~ t& in this way. We will also be 
able to use observed values of p and p\ to constrain the two-dimensional multiverse 
forces on these parameters. We will find, however, that these empirical constraints 
lead to tension with theoretical expectations. This, therefore, disfavors the measure we 
consider here. 

In the next section we will consider a different class of measures which will be 
more successful, and also more powerful in that they constrain an additional scanning 
parameter. 

5.1 Catastrophic boundary from large scale structure 

The role of the measure In theories with at least one long-lived metastable vacuum 
with positive cosmological constant, a finite false vacuum region can evolve into an 
infinite four-volume containing infinitely many bubbles of other vacua. Each bubble 
is an infinite open universe; if it contains any observers, no matter how sparse, it will 
contain infinitely many. These infinities must be regulated if we are to say that some 
outcomes of observations are typical and others unlikely. This is the measure problem 
of eternal inflation. 

We will consider four measure proposals that are reasonably well-defined and do 
not suffer from known catastrophic problems. 6 The first and second are the causal 
patch [5] (CP) and causal diamond (CD) [6] cutoffs; the third and fourth are versions 

5 For definiteness, we will consider only positive values of the cosmological constant, although it is 
not difficult to extend our analysis to negative values. 

6 The measures used in the older literature on environmental selection of the cosmological constant, 
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of the scale factor cutoff, which we will refer to as SF1 [8] and SF2 [9]. We will not 
review these measures here in detail, but focus instead on their implications for the 
analysis of the cosmological constant problem. The measures differ in important ways: 
they lead not only to different catastrophic boundaries, but also to different multiverse 
forces on the cosmological constant and virial density. 

Both measures are based on a geometric cutoff: one computes the expected number 
of observations made in some finite spacetime region. They differ in how that region 
is defined. The causal patch measure (CP) restricts to the causal past of the future 
endpoint of a geodesic in the spacetime. In long-lived metastable vacua with posi- 
tive cosmological constant, this definition reduces to counting observations that take 
place within the cosmological (de Sitter) event horizon surrounding the geodesic. This 
has two important consequences. First, observations made after vacuum domination 
(^obs > *a) contribute negligibly; they are exponentially suppressed by the accelerated 
expansion that empties the horizon region of matter. This leads to a catastrophic 
boundary A < t~^ s . 7 Secondly, at all times before vacuum domination (t < £a), the 
total matter mass inside the causal patch scales as A -1 / 2 . For smaller A, more matter 
is included, which mitigates the landscape force towards large values of A. 

The causal diamond (CD) cutoff further restricts the surviving spacetime region. 
A point is included only if it lies both in the causal patch and in the causal future of 
the intersection of the generating geodesic with the reheating hypersurface. This leads 
to the same catastrophic boundary, A < i~ b 2 s . However, before t\, the comoving volume 
of the diamond is controlled by the future light-cone. Thus it is independent of A and 
grows linearly with time. This favors later observers. 

As shown in Ref. [9], the scale factor cutoff (SF1) can be defined as a small neigh- 
borhood, of fixed physical volume, of a geodesic in the multiverse. (The formulation of 

such as observers-per-baryon, fail catastrophically in eternally inflating universes driven by false vacua 
[17,18], so they will not be considered here. For other examples of catastrophic problems, see, e.g., 
Ref. [19]. 

7 This boundary is stronger than the catastrophic boundary coming from the disruption of galaxy 
formation. Of course, nothing prevents the formation of observers in a universe where suitable galaxies 
form. However, such observers are not counted if the galaxies are not contained in the cutoff region. 
The existence of such a boundary, coming entirely from the measure, may seem counter-intuitive; in- 
deed, it is among the most interesting lessons learned from recent studies of the measure problem [6] . 
In eternal inflation, no conventional anthropic boundary is well-defined without a measure. For ex- 
ample, with A 3> t~ir> there is still a tiny amplitude for the formation of galaxies from unusually 
strong density perturbations, so there will be infinitely many observers no matter whether A is above 
or below the Weinberg bound. To make a meaningful comparison, we need to pick a measure. The 
measures usually considered are geometric cutoffs that necessarily depend on quantities that affect the 
geometry, such as A. It should not surprise us, then, that they can introduce catastrophic boundaries 
on observable parameters. 
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the scale factor measure given in Ref. [8] is not completely equivalent but shares the 
following properties [9] 8 relevant for our analysis.) Like in the CD measure, and unlike 
CP, the size of the SF1 cutoff region is independent of A before the time t\. Because 
the geodesic is likely to become trapped in a collapsed region, and one is including only 
its immediate neighborhood, the expansion of the universe after £ vir does not reduce 
the number of observers contained in the cutoff region. This leads to two important 
differences from both CD and CP. First, there is no suppression for t Q b s > t\, so the 
only relevant catastrophic boundary is the Weinberg bound, A < t^?. And secondly, in 
SF1 the number of observers is proportional to their physical density near the geodesic; 
this favors larger virial density. 

SF2 is a hybrid of SF1 and CD. This measure has only been tentatively formu- 
lated [9], but we include it here because we are mainly interested in the consequences 
of a measure with the following properties: like in the CD measure, the catastrophic 
boundary is A < t~ b 2 s , and the spatial size of the cutoff region is A-independent before 
t\. The number of observers is proportional to the average density of matter in the 
universe at the time when observations are made. This is similar to SF1, except that 
for SF1 the relevant density is set at an earlier time, when the first collapsed structures 
form. 

In this section, we begin by considering the catastrophic boundary relevant to SF1: 
the disruption of the formation of galactic halos. In section |6], we will consider the 
catastrophic boundary relevant to CP, CD, and SF2: the exponential dilution of the 
number density of galaxies. For each case, we will obtain predictions for scanning 
parameters and empirical constraints on multiverse forces. We will discuss whether 
these constraints are compatible with theoretical expectations, which also depend on 
the measure. 

The boundary In a given universe, the cosmological constant disrupts the formation 
of halos above a certain mass, because the logarithmic factor f(M) in Eq. ( |3.6|) will 
enter into the condition of Eq. ( |3.15|) . It is possible that there is a minimum galaxy mass 
necessary for the formation of observers. As in the previous section, however, we will 
work with the conservative assumption that any galaxies will do. Thus, a catastrophic 
boundary for the cosmological constant is reached only when is so large as to disrupt 
the formation of the smallest halos that can cool. By Eq. ( |3.15[ ), this corresponds to 
the condition 

p vir (M_(p)) 

PA,max(p) = , (5.1) 

54 



8 This was not originally recognized in Ref. [8], who concluded that the properties of their formu- 
lation were similar to those ascribed here to SF2. 
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Figure 5: When both p\ and the virial density p scan, the resulting behavior depends on 
the direction of the multiverse force and on the boundaries. On the left, a multiverse force 
with only a p\ component leads to runaway behavior along the virialization boundary, as 
shown by the dashed arrow. On the right, a component of the force in the direction of p 
is added, as well as the boundary from galactic radiative cooling. Typical universes lie at 
the tip of the cone, indicated by the dot. Strong forces would locate the dot within a factor 
of two of each boundary. But the forces may be mild, placing the dot one to two orders of 
magnitude from the boundaries, as observed in our universe. With an accuracy that depends 
on the strength of the force, multiverse predictions result for both p\ and p in terms of the 
fundamental parameters that determine the location of the boundaries. 



where p vir (M_) and M_(p) are given in Eqs. ( pT6| ) and ( |3.10| ). 

Up to logarithmic corrections from /(M_), this catastrophic boundary corresponds 
to a line of slope 1 in the (p, p\) plane shown in Fig. |5]. Of course, the analysis in the 
following subsection would apply to any measure in which this line is the relevant 
catastrophic boundary. 

5.2 Predictions 

We have already assumed a force towards large vacuum energy p\, and we have estab- 



lished two catastrophic boundaries, Eqs. (|4.1| , |5.1|) . These boundaries intersect at the 
critical point 

(p,Pa) = (PoPa,c), (5.2) 



where 



PA, C = PA, max (p c ) = ( %yl}f ) - (2-8 x ID" 3 eV) 4 . (5.3) 



Here, we have used the expression for p vir c given in Eq. ( |4.3| ). 
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Whether these critical values become predictions depends on the sign and magni- 
tude of the force component in the p direction, pp = din f /din p. (Note that this is 
not the same as the force on p considered in the last section, which was understood to 
result from integrating out all other scanning parameters, such as A.) 

The left panel of Fig. [| shows the runaway problem that would arise in the absence 
of a sufficiently strong force toward small p, i.e., if pp > — Pa- The total force points to 
the structure formation boundary, but once this boundary is reached, the component 
tangential to the boundary goes in the positive pa and p directions, giving the runaway 
behavior indicated by the dashed line. The cosmological constant problem is not solved 
in this case. Moreover, the prediction for p in the previous section would also be lost, 
since the effective force on p, after integrating out A, would point towards large values. 

Let us assume, therefore, that pp < —p\, i-e., the force towards small p is strong 
enough to overwhelm the force towards large pa once the catastrophic structure forma- 
tion boundary is reached. In this shown in the right panel of Fig. [|, the force 
selects the tip of the cone formed by the two boundaries as the most likely universe, 
and we obtain two predictions: 

p ~ p c ~ (0.97 x 1(T 4 eV) 4 , 

(5.4) 

PA ~ PA,c ^ (2.8 x 1(T 3 eV) 4 . 

Recall that p c and p v ir,c, given in Eqs. ( |4.1[ ) and (|4.3|), have the same parametric 
dependence on fundamental parameters and differ only through numerical factors. Ne- 
glecting such factors, we are able to relate both p and pa to fundamental parameters 
as 

G N m 4 e m 2 

P~Pa~ j • ( 5 - 5 ) 



As a corollary, we are able to explain the double coincidence 

ft 2 

tA ~ trad ~ ^vir ~ ^ • (5.6) 

Thus, we predict that most observers will find themselves in a universe with p 
larger, but not much larger, than p c , and with pa comparable to p v i r ,c/54. (Note that 
Pa could be smaller or larger than this critical value, because the catastrophic boundary 
is not a line of constant pa- However, it should not differ very much from the critical 
value.) 

The first prediction in Eq. (|5.4j ) reproduces our successful result in section fy| and 
it yields the same corollaries: predictions for the characteristic mass, virial density, and 
cooling timescale of galactic halos. The second prediction addresses the cosmological 
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constant problem. This prediction, too, is successful: the observed vacuum energy, 
Pa.o — (2.3 x eV) 4 , is only a factor 2.1 smaller than the critical value. 

This might suggest that the force on the cosmological constant must be strong. But 
this is not the case. To gauge the multiverse forces, we must compute the distances 
orthogonal to the catastrophic boundaries. We shall now see that this constrains the 
force on the cosmological constant to be quite weak. 

5.3 The strength of the multiverse force 

We now have a two-dimensional parameter space, with variables 

(5.7) 
(5.8) 





= In 


P_ 






Pc 




= In 


PA 






PA,c 



By Eq. (|2.7|), the matrix 



defines two new parameters 



ui = In 4- , (5.10) 

Pc 

, PA, max IP J /c i 1 ^ 

u 2 — In — ■ , (5-11) 

PA 

that are orthogonal to catastrophic boundaries at U\ — and m 2 = 0, where pA, m ax is 
given in Eq. fl5.1|). (We have used that pA,max is proportional to p, up to a logarithmic 
term which we neglect.) 

By Eq. ( |2.9| ), (i? _1 ) T defines associated force components S, with 

si = Pp+PA, (5.12) 
s 2 = ~Pa, (5.13) 

acting on the parameters U . The absence of a runaway problem corresponds to the 
condition si < 0, i = 1, 2, which agrees with the force conditions we identified in 



section 5.2 above 



The observed values in our universe are 



Wi jO ~ln40 -3.7, (5.14) 
u 2 , ^ In 139 ~ 4.9. (5.15) 
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The second line is the statement that the cosmological constant in our universe is more 
than a hundred times smaller than necessary for large scale structure formation. By 
demanding that these observed values are within one or two standard deviations from 
the median, we obtain from Eq. ( |2.13| ) that 



P? + PA = -[0.19lg5i(la)lSS(2(T)], (5.16) 
-p A = - [0.141 tS(la) i£gg(2*)] . (5.17) 

Note that the right hand side of Eq. ( |5.16| ) is the same as the constraint on the force 
Pp derived in the one-dimensional setting of Eq. ( f4.9|) . This is as is should be. The left 
hand side looks different, but this is a notational problem: in the previous section, pp 
was understood to include all contributions from integrating out other parameters and 
boundaries. Here, however, the contribution from p\ was excluded in the definition of 
Pp. It thus appears explicitly in Eq. Q5.16|) . 



How do these bounds compare to theoretical expectations? Because is not a 
special point in the spectrum of the cosmological constant, one can treat the statistical 
distribution of p\ as flat in a small neighborhood of (and only an extremely small 
neighborhood is of interest anthropically) : 

■ const, (5.18) 



dp A 

at least after smoothing over intervals much smaller than the observed value. Thus, 
the simplest theoretical expectation is a force of unit strength towards large values of 
A: 

PA = 77"j = I- 5.19 

o In p A 

This result, however, is modified (in different ways) by all of the measures we consider 
here. 

The catastrophic boundary we have used — the Weinberg bound on the formation 
of galactic halos — arises from the scale factor measure as formulated by De Simone et 
al. [8,9]. Correspondingly, we should here consider the implications of this particular 
measure for the strength of F\. We will compare this theoretical expectation to our 
empirical constraint, Eq. ( p. 17] ). 

Let n D bs be the number density of observers per physical volume at the time t v ir+^obs 
when observations are made in a given bubble universe. As shown in Ref. [9], the 
scale factor measure weights observers by the physical number density these observers 
would have had at the time when the first structures formed that would eventually be 
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incorporated into their halo. This weight is approximately 

w SF1 oc ( J p vir , (5.20) 

VPobs/ 

where p Q b s is the matter density at the time t v - ir + t ^ s . The first factor is the number of 
observers per unit mass. In our sharp-boundary approximation, this factor is constant 
and nonzero on the anthropic side of the catastrophic boundaries in the (p, pa) plane. 
Beyond the cooling and structure boundaries, it vanishes. 

The second factor, up to numerical coefficients, agrees with p, so the weight wsfi 
contributes a single power of p to the distribution function /. Thus, the force on the 
anthropic side receives a contribution 

(Ap p -,Ap A ) SF1 = (l,0) (5.21) 

from the measure. Therefore, the total force is (pp,p\) = (1 + Pp, 1), where pp = 
din f /din p is the unknown contribution from the distribution of p among landscape 
vacua. For the linear combinations constrained in Eqs. (|5.16 , |5.17| ), we find 



(Pp + Pa,Pa) = (2+P p -,1). (5.22) 

The first constraint can be satisfied by assuming that the unknown prior landscape 
force, pp, lies between —2 and —3. While this may seem like an implausibly strong force 
towards small virial density — an inverse power law of mass dimension —8 to —12 — we 
cannot rule out that it may be realized in the string landscape. 

The second constraint, however, on the force on A, is definitely incompatible with 
the theoretical expectation at the 2a level. (It would be compatible at the 3a level.) 
Thus, with the scale factor measure in the formulation of De Simone et al. [8] , multiverse 
forces cannot simultaneously explain p vir and p\ without additional assumptions. 



6. Predicting the Cosmological Constant and Observer Time 
Scale 

In this section, we will again consider the case of scanning both the cosmological con- 
stant pa and the virial density parameter p, but using a different measure for regulating 
the divergences of the eternally inflating multiverse. In the previous section we used 
the scale factor measure in the formulation of Ref. [8] (SF1). We will now explore the 
implications of the causal patch (CP) [5], causal diamond (CD) [6], and modified scale 
factor (SF2) [9] measures. We will find important differences to the previous section, 
and subtle differences among the three measures we now consider: 
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• All three measures yield the same set of three catastrophic boundaries, in the 
three-dimensional parameter space (p, PA,t bs)- The first two boundaries are the 
same as for the SF1 measure; the boundary on t obs is new. 

• Because of the appearance of t Q ^ s , we will allow three parameters to vary in this 
section: p, p\, and t G b s - 

• The three measures make different contributions to the multiverse force vector. 
We will find observational constraints on the corresponding multiverse forces that 
are in excellent (CP, CD) or reasonable (SF2) agreement with theoretical expec- 
tations. 

6.1 Catastrophic boundary from galaxy dispersion 

All three measures considered in this section exponentially suppress the weight of ob- 
servations that take place after the cosmological constant dominates, like exp (— 3t/tj\). 
There is some ambiguity about which critical value to pick in the sharp-boundary 
approximation. We will use t A = (87rG N p A /3)- 1/2 , but our results do not depend 
sensitively on small modifications of this choice. 

Thus, in the sharp boundary approximation, we require 

tvir + ^obs < ^a, (6.1) 

as a necessary condition for the existence of observers. Here, we define t b s as the 
time delay between the formation of a galactic halo at time t v i r and the time when 
the halo hosts observers. (We neglect the cooling time scale, which must not exceed 
t v i r in any case.) The delay t Q ^ s includes the time needed for the formation of a solar 
system and Darwinian evolution, or whatever other forms the development of life in 
gravitationally collapsed structure may take. We have no idea how to relate t Q ^ s to 
fundamental particle physics or cosmological parameters, but we will assume that it 
depends on sufficiently many parameters, sufficiently sensitively, that we can treat it 
as scanning independently of the other four parameters we are varying in this paper 
(p, t r ad, Pa, and Q). 9 

In the presence of power law forces on p and t b s , we make a negligible error if we 
write Eq. (|6.1| ) as two independent catastrophic boundaries: 

tobs < ^a, (6.2) 
t vir < t A . (6.3) 

9 Note that our treatment of t bs is not very different from t\. Neither quantity can be simply 
computed from "fundamental" parameters of our vacuum, but both are expected to scan in the 
multiverse. We predict both quantities by correlating them with other scales. 
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The existence of galaxies, of course, remains a necessary condition for observers, and 
it yields a bound tighter than Eq. ( |6.3| ): from Eq. ( |5.1|) we find 

ivir < 7^t A ~ 0.60t A . (6.4) 

Thus, we will use Eq. fl5.2|) as the new catastrophic boundary for the new parameter 
Wi- We will continue to treat p v ; r /54, from Eq. ( |5.1|) , as a catastrophic upper bound 
on pa, and p c , given in Eq. (|4.1|) , as a catastrophic lower bound on p. 



6.2 Predictions 

The three boundaries we have established intersect at the critical point 

(p, Pa, t ohs ) = (p c , pa, c , tobs.c), (6.5) 

where 

W = U(pa,c) = 1/2 Q 2 - 11-0 Gyr. (6.6) 

3 7 T C ion Aad ^^ e m p 

Supposing that the three-dimensional multiverse force points towards this critical 
point, we predict that observed values of these parameters should be within one or two 
orders of magnitude of the critical values. The predictions for p and pa are the same as 
in the previous sections. The new prediction for t Q ^ s is that it should not differ much 
from t bs,c, given above: 

tabs ~ £obs,c- (6.7) 

The observed value is 

t obs , D = 10.1 Gyr, (6.8) 

so this third basic prediction is also successful. 
We have now explained the triple coincidence 

a 2 

tabs ~ *A ~ ^rad ~ ^vir ~ ~^ 5 • (6.9) 

Next, we will determine under which conditions the multiverse force points towards 
the critical point, and we will establish empirical constraints on the the three force 
components. 

6.3 Forces with the causal patch measure 



In addition to the variables z\ and z 2 defined in section |5.3| , we now have the third 
scanning parameter 

z 3 = ln-^. (6.10) 

£obs,c 
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The matrix 

/l 

5=1-1 J (6.11) 
\0 -1 -2 



defines the boundary-orthogonal parameter 

3 1 

8ttGn PAi obs 

By Eq. fl2.9|), (-B _1 ) T defines a boundary-orthogonal force vector S, with 



"3 = In ( -it- J • ( 6 - 12 ) 



si = Pp + Pk - ^Pobs, (6.13) 
1 / 

52 = "PA + ^ Pobs, (6.14) 

1 / 

53 = -^Pobs, (6.15) 

where p Q b s = dln//01ni o b s - The absence of a runaway problem, i.e., the condition that 
the multiverse force points towards the critical point, is that all < 0. In particular, 
we learn that it is necessary to assume a statistical preference for a large time scale for 
the development of observers in gravitationally bound regions. This assumption seems 
plausible, given the complex processes that are likely to be involved. Since p\ appears 
only in the combination p\ — |p b s , this preference acts counter to the multiverse force 
towards large cosmological constant. This, too, is intuitively plausible: there are fewer 
vacua with smaller cosmological constant, but on average they contain more observers. 
From the observed values in our universe we find 

Ui )O ~ln40 ~3.7, (6.16) 
u 2i0 ^ In 139 ~ 4.9, (6.17) 
m 3i0 ~ In 2.5 ~ 0.92. (6.18) 

The fact that u% is of order unity is closely related to the coincidence problem: we 
live in the unique era when vacuum and matter energy densities are comparable. By 
demanding that these observed values are within one or two standard deviations from 
the median, we obtain from Eq. (|2.13|) that 

pp + pa - ^Pobs = - [o.i9 ±8:?l(io-) ±8:f|(2o-)] , (6.19) 

~Pa + ~Pob s = - [0.141 jKjj(la) ±° itt(2a)] , (6.20) 
— Pots = - [0.76 + _lf 7 (la) + J 3 7 l(2a)] . (6.21) 
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The first two constraints are identical to those found in the previous section, except 
that the contribution to the multiverse force from p Q i, s is now made explicit on the left 
hand side. 

How do these bounds compare to theoretical expectations? As reviewed in the 
previous section, the statistical distribution of A in the landscape is flat, giving a prior 
force Pa = 1 towards large values. As discussed at the beginning of this section, the 
three measures differ in the spatial size of the cutoff region. On the anthropic side of 
the catastrophic boundary, the number of observers included depends differently on t b s 
and pa for each measure. This means that the measures make different contributions 
to the multiverse force vector. We will discuss each measure in turn, beginning with 
the CP measure. 

The CP measure contributes a factor p~^ 2 [6]: 

^cp = (^) Pl l/ \ (6.22) 

V Pobs / 

where the first factor, the density of observers per unit mass, is again approximated as 
constant on the anthropic side of the catastrophic boundaries and zero beyond. Thus, 
the CP measure makes no contribution to the forces on p and t b s , and we find: 

(Ap- P , Ap A , Ap ohs ) CP = (0, ~, 0). (6.23) 

The CP measure, therefore, leads to a reduced force on the cosmological constant: 

PA = \ (6-24) 

In Eq. ( |6.20| ), p\ and the force towards large observer time scale appear with opposite 
sign, reducing the magnitude of the left hand side further. This should be compared to 
the analogous equation Eq. ( |5.17j ) for the SF1 measure, which contained only p A = 1, 
and so could not be satisfied at the 2a level. By contrast, the above three empirical 
constraints are all compatible with Eq. (|6.24j ) at the la level. Possible values for 
(Pp,Pobs) include (-§, §) and (-±, §). 

6.4 Forces with the causal diamond measure 

How would this discussion be modified if we used the CD measure? The catastrophic 
boundaries remain unchanged, preserving both the general prediction that we should 
find ourselves near the triple critical point (p c , p\ tC , t bs,c), and the empirical constraints 
given in Eqs. (|6U9| - |OlD . 



However, the measure contributes differently to the multiverse force. It does not 
contribute to the force on p\, but instead favors late times of observation, t vir + t ohs , 
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because the mass inside the causal diamond depends linearly on this quantity as long 
as it does not exceed t\. This corresponds to a contribution of a force of strength 1 on 
the variable t vir + t G b s : 

Ap vir+ obs = 1- (6.25) 

In terms of the individual variables, t v ; r and £ b s , this implies that for t va ^> t bs, the CD 
measure contributes (App, Ap^ Ap obs ) CD = (— |, 0, 0); in the opposite limit, (0, 0, 1). 

Thus, the force is not uniform in the allowed quadrant of parameter space; instead, 
it changes across the hypersurface t v - ir ~ t G bs- This case cannot be treated with the 
simple analysis following Eq. ( |2.8|) , and we will not attempt here to constrain the 
prior force strengths and directions compatible with the observed values of the three 
parameters. 

We note, instead, that the measure force Ap vir+0 b s acts to increase t v i r + t Q bs- Be- 
cause of the constraint Eq. Q6.1J ) that arises in this measure, it also acts to increase t A 
with strength 1, or equivalently to decrease A with strength —1/2. Therefore, the mea- 
sure force effectively counteracts the prior pressure towards large A, reducing it from 1 
to 1/2. This is similar to the causal patch, where this reduction happens more directly, 
without the need to invoke one of the catastrophic boundaries; see Eqs. ( |6.23| , |6.24[ ). 
We thus expect that the causal diamond, like the patch, leads to a good agreement 
between observed and predicted values for these three parameters, for a reasonable 
range of prior force strengths (j>p,p b s )- 

6.5 Forces with the modified scale factor measure 

Finally, let us turn to SF2, a version of the scale factor measure that requires t v - ir + 
£obs < t\ [9]. We may again translate this condition into the two boundaries given 
in Eqs. ( |5.1| , |6.2|) . 10 Thus, the catastrophic boundaries are the same as for the other 
two measures considered in this section. And so the SF2 measure, too, gives rise 
to the general prediction that we should find ourselves near the triple critical point 
(p c , Pa,c, ^obsx), and to the empirical constraints given in Eqs. ( |6.19| - |S~2"T|) . 



However, the SF2 measure makes contributions to the forces that differ strongly 
from those arising from the CP or CD. It weights by the physical density of ob- 
servers, which is proportional to (t v ; r + t bs)~ 2 - Thus, for t vir t Q b s , it contributes 



10 Strictly, the SF2 measure involves an effective i b s related to the time when the latest geodesic 
reaching observers achieved maximum expansion. One would expect that this occurs at about one 
half of the time t v ; r + i b s , so one should redefine i bs,SF2 = (^obs.cp — t v i r )/2. However, compared to 
the CP, the SF2 measure suppresses large values of pa > ^obs/^ N an additional factor p A 3 ^ 2 ; so 
the sharp boundary on t bs, which was iobs.cp < t\, should also be moved to a smaller value. If we 
choose £ bs,SF2 < £a/2, then these two effects roughly cancel, so we will omit both from our analysis 
for simplicity. 
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(App, ApA, Ap bs)sF2 = (1) 0, 0); in the opposite limit, (0, 0, —2). This is unhelpful, be- 
cause both contributions go in the wrong direction: by Eq. ( |6.21|) , p D b s must be positive, 
and by Eqs. ( |6.19| , |6.20| ), pp must be negative. Thus, the empirical constraints require 
that large landscape forces pp < and p D b s > overcompensate for the contribution 
from the measure. This is particularly acute for p ^ s , which must be larger than in the 
CP measure by about 1, in order to compensate for the larger value of p\ (1 instead 
of 1/2) in Eq. 06.2OQ . But we cannot exclude the SF2 measure; it is possible that the 
prior distributions in the landscape satisfy these constraints. 

7. Predicting the Primordial Density Contrast 

Our halo, which is among the largest galactic halos, virialized approximately 3 Gyr 
after the big bang. According to Eq. ( |3.21|) , Compton cooling in our universe became 
ineffective only a little earlier, at 360 million years. We have argued that there is 
an environmental understanding of why the Milky Way halo was just able to cool 
radiatively, could it be that there is also an environmental reason to explain why it just 
failed to cool by inverse Compton scattering? 

Indeed, the proximity to Compton cooling is striking. For example, with Q a factor 
of 4 larger than in our universe (and holding fixed the temperature at equality), a halo 
with the mass of the Milky Way would have undergone Compton cooling. This suggests 
that there is a catastrophic boundary associated with Compton cooling, and that the 
multiverse distribution has pushed us close to it. In this section, we will explore the 
consequences of such a boundary. We stress that we do not explain why the anthropic 
weighting factor varies suddenly across the Compton boundary; we will simply pursue 
the consequences of assuming that it does. 

7.1 Catastrophic boundary from Compton cooling 

As discussed in section ^74], Compton cooling is effective if t com p < £ v irj this is the case 
for all galaxies that form prior to the cosmological time t C omp,max given in Eq. (|3.21 ). 
We are assuming that Compton cooling inhibits the formation of galaxies, or changes 
their structure sufficiently to make them inhabitable. In our universe, Compton cooling 
only affects halos that virialize very early and thus have small mass. We have defined a 
catastrophe as the absence of habitable galaxies of any mass scale. Thus, a catastrophe 
occurs if parameters are altered so that Compton cooling affects all halos, including 
those of the largest mass able to cool radiatively, M + , given in Eq. ( |3.14| ). Namely, a 
catastrophe is averted if 

^comp,max *^ ^vir (-^+) ■ C^ - l) 
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Using Eq. (|3.4| ), we thus find the catastrophic boundary 



2 10 7t 1 /3/(m + ) 5 a 4 

We could have expressed this boundary in terms of any two of the three parameters 
Q, T eq , and p = Q 3 T* . For continuity with the previous sections, we have chosen to 
retain p as a scanning parameter. 

There is a subtlety here. Compton cooling can be catastrophic only if it is faster 
than radiative cooling; otherwise, it would have no opportunity to act. However, this 
condition will automatically be satisfied if the above inequality is violated. For the 
largest halos able to cool radiatively, t rad = £ vir by definition, so a violation of Eq. ( [7.1[) 
implies that Compton cooling is faster than radiative cooling. For smaller halos, which 
virialize earlier, t comp oc t^ 3 whereas t rac \ oc t\[ 3 . Therefore, Compton cooling will be 



faster than radiative cooling for all halos capable of cooling if Eq. (|7.1]) is violated. 
7.2 Predictions 

We can now consider the four- dimensional parameter space (p, p\, t b s , Q)- We have 
established four catastrophic boundaries, which intersect at the critical point 

(p~c> PA,c, ^obs,c, Qc)- (7.3) 

The first three critical values have been derived in the three previous sections; the 
fourth can be obtained by setting p — > p c in Eq. ( [7.2|) . Using Eq. ( |4.1| ), we find 

„ 7r 1//3 c v i r ^ co iA^ eQ f D Kfn m e q 

Q c = Z J — ~ 2.4 x llT 3 . 7.4 

As usual, let us now assume that the four-dimensional multiverse force is directed 
towards the critical point. (We will analyze this condition in detail in the next subsec- 
tion.) Then we predict 

Q ~ Qc (7.5) 

in addition to the three other predictions (p,p\,t ohs ) ~ (p c , p\ c , t ohsc ) established in 
the previous sections. The observed value is Q a ~ 2.0 x 1CT 5 , a factor of 120 below the 
critical value. As a corollary, we are now able to explain the quadruple coincidence 

a 2 

£comp,max ~ ^obs ~ ~ ^rad ~ *vir ~ 5 (7.6) 

of timescales observed in our universe. 
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7.3 Forces 



In addition to the variables Zi,Z2,z% denned in section |6.3| , we now have the forth 
scanning parameter 



The matrix 



B 



( 1 \ 
1-10 
0-1-2 
\-\ -2 J 



defines the boundary-orthogonal parameter 



(7.7) 



(7- 



114 = In 



QlPc 
Q 2 p 



(7.9) 



and leaves the definitions of Ui,u 2 , u 3 unchanged compared to section p73 . 
By Eq. ( |2.9|) , (5 _1 ) T defines a boundary-orthogonal force vector S, with 



si = Pp + Pa ~ 2 P° hs ~ 2 PQ ' 



S-2 
S 3 



'PA + -Pobs, 
1 

2 
1 



" 2 Pobs 5 



S 4 = --PQ- 



(7.10) 

(7.11) 
(7.12) 
(7.13) 



The absence of a runaway problem, i.e., the condition that the multiverse force points 
towards the critical point, is that all s, < 0. In particular, we must assume a statistical 
preference for large values of Q, which seems entirely natural. 



(7.14) 
(7.15) 
(7.16) 
(7.17) 



our 


universe 


we find 






In 40 


~ 3.7, 


U2,o 




In 139 


~ 4.9, 


«3,o 




In 2.5 


~ 0.92, 


U4,o 




In 355 


~ 5.9. 
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By demanding that these observed values are within one or two standard deviations 
from the median, we obtain from Eq. (|2.13|) that 

Pp+Pa- ipob. -\pq = - [0.19 t° H(la) t° f 8 (2a)] , (7.18) 

-PA + iflob. = " [0.141 ™(l<x) 1S(2,)] , (7.19) 

-\poU = - [0.76 tl§(la) iS$(2*)] , (7.20) 

= - [0-12 l°;*j(la) 18;g(2a)] . (7.21) 

The first three constraints are identical to those found in section |6.3| , except that the 
contribution to the multiverse force from pq is now made explicit on the left hand side. 

How do these bounds compare to theoretical expectations? We will discuss only 
the causal patch measure here, which leads to the force p\ = | on the cosmological 
constant. With this value, the above constraints can easily be satisfied at the la 
level with natural force strengths, e.g., with (pp,p bs,PQ) = (— f, §, |), (— f, §, |) or 

V 4? 4' 2/" 



8. Discussion 



We have found that multiverse forces and catastrophic boundaries can explain the 
quadruple coincidence of timescales, Eq. (p..l|), or equivalently, the observed values of 
the virial density, the cosmological constant, the observer timescale, and the timescales 
for radiative and Compton cooling. 

As more parameters are taken to scan, dare we think that all fundamental param- 
eters could be determined in this way? The answer is clearly no: as far as we know 
many parameters, such as the strong CP parameter and the r lepton mass, are not 
relevant for any catastrophic boundary. Still, could it be that all the parameters that 
do enter catastrophic boundaries get determined environmentally? For the boundaries 
considered so far, the relevant set of parameters is 

a,m e ,m p ,G N ,Q,T eq ,p A ,r] B , (8.1) 

where T)b is the number density of baryons relative to photons. Depending on the 
measure used to compute the statistical averages we could also include t b s - If all the 
parameters of Eq. ( |S.1D scan independently then the structure and cooling boundaries 
can determine p\, Q, T eq (and t ohs ), but what of the others? 

The ratio m e /m p plays an important role in the boundaries of nuclear stability, 
leaving a,m p ,r]B and Gn- The proton mass is directly derived from the QCD scale, 
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Aqcd, and we can assume that this is related to a by gauge coupling unification 11 — i.e. 
we are not quite right to say that all the parameters of Eq. (|8.1|) scan independently. 
Taking to set the unit of energy, Mpi, leaves just two undetermined dimensionless 
parameters: t)b and Aqcd/-^pi- Thus the landscape force and catastrophic boundaries 
have rigidly determined the entire parameter set of Eq. ( |8.1| ), with the exception of the 
baryon asymmetry and an overall mass hierarchy. In fact, they may overdetermine it 
if t bs is a function of the parameters listed in Eq. flS.l|) . This, however, is not the case 
if t hs is very sensitive to the values of a, m e , ■ • • , so that it can take a wide range of 
values within the parameter space in which the fundamental parameters, a, m e , • • • , 
satisfy the nuclear stability and other requirements. In such a case, t Q b s can be viewed 
as an independent scanning parameter. 

Our universe lies at the tip of a cone in the four dimensional parameter space of 
{m e /m p , Q, T eq , p^} (or of five dimensional space including £ bs); but when an additional 
dimension such as tjb and Aq CD /M P i is added, a runaway behavior could reemerge. It is, 
however, possible that the direction opened by the addition of a new parameter is again 
constrained by a new catastrophic boundary containing that parameter. For example, 
the up and down quark masses are constrained by further nuclear stability boundaries, 
and the weak scale and baryon asymmetry play a crucial role at the boundary associated 
with big bang nucleosynthesis. Given these possible boundaries, we can imagine that 
most of the environmentally relevant parameters in the standard particle physics and 
astrophysics are determined by multiverse physics as described here. 

Ultimately, an important question will be what determines the overall mass hier- 
archy between the Planck and particle physics mass scales, allowing the existence of 
large observable universe. We see three different possibilities: 

• The multiverse distribution function. For example, the unified coupling at the 
Planck scale may have a distribution f(a) peaked around some perturbative value 
that corresponds to a peak in the distribution for Aqcd exponentially below Mp\. 

• A further catastrophic boundary. As Aqcd/-^pi grows, complex structures such 
as the sun will get smaller, and eventually cause catastrophic changes. Perhaps 
our universe is already close to some such boundary. In this case, n catastrophic 
boundaries constrain n scanning parameters. (Throughout this paper we have 
assumed a multiverse force that places typical observers in universes near the 
tip of a cone. An alternative possibility is that the observer region in the n 
dimensional space is very small.) 

n The relation involves the weak scale, which may be determined by some other boundary. 
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• The discretuum. The vacuum energy scales as the 6th power of the overall hi- 
erarchy: pa oc (Aq CD /M p1 ) 6 , given that environmental selection by nuclear sta- 
bility sets m e ~ m p ~ Aq CD . Suppose that f(a) gives a force to low values 
of (Aqcd/^pi)- The vacuum energy may become so small that the discreteness 
of the string landscape will become apparent, eventually stopping the potential 
runaway. If this is the origin of the overall hierarchy, the number of string vacua 
is of order M^/p A « 10 120 . 

At present we are unable to compute the overall hierarchy, no matter how it is de- 
termined. Nevertheless, we are encouraged that simple arguments based on galaxy 
cooling allow a broad numerical understanding of p\, M gal , Q and T eq in terms of the 
fundamental parameters in particle physics, a, m e , m p and G^. 
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